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Abstract. We study the geometric phase of an open two-level quantum system under the influence of 
a squeezed, thermal environment for both non-dissipative as well as dissipative system-environment in- 
teractions. In the non-dissipative case, squeezing is found to have a similar influence as temperature, of 
suppressing geometric phase, while in the dissipative case, squeezing tends to counteract the suppressive 
influence of temperature in certain regimes. Thus, an interesting feature that emerges from our work is the 
contrast in the interplay between squeezing and thermal effects in non-dissipative and dissipative interac- 
tions. This can be useful for the practical implementation of geometric quantum information processing. 
By interpreting the open quantum effects as noisy channels, we make the connection between geometric 
phase and quantum noise processes familiar from quantum information theory. 



PACS. 03.65.Vf Phases: geometric; dynamic or topological 
03.67.Lx Quantum computation 



03.65.Yz Decoherence; open systems 



1 Introduction 

Geometric Phase (GP) brings about an interesting and im- 
portant connection between phase and the intrinsic cur- 
vature of the underlying Hilbert space. In the classical 
context it was introduced by Pancharatnam [T], who de- 
fined a phase characterizing the intereference of classical 
light in distinct states of polarization. Its quantum coun- 
terpart was discovered by Berry [2] for the case of cyclic 
adiabatic evolution. Simon |3] showed this to be a conse- 
quence of the holonomy in a line bundle over parameter 
space thus establishing the geometric nature of the phase. 
Generalization of Berry's work to non-adibatic evolution 
was carried out by Aharonov and Anandan and to the 
case of non-cyclic evolution by Samuel and Bhandari [5] , 
who by extending Pancharatnam's ideas for the interfer- 
ence of polarized light to quantum mechanics were able 
to make a comparison of the phase between any two non- 
orthogonal vectors in the Hilbert space. An important de- 
velopment was carried out by Mukunda and Simon [5], 
who, making use of the fact that GP is a consequence of 
quantum kinematics, and is thus independent of the de- 
tailed nature of the dynamics in state space, formulated a 
quantum kinematic version of GP. 

Uhlmann was the first to extend GP to the case of 
non-unitary evolution of mixed states, employing the stan- 
dard purification of mixed states T|. Sjoqvist et al. [8] 
introduced an alternate definition of geometric phase for 
nondegenerate density opertors undergoing unitary evolu- 
tion, which was extended by Singh et al. W to the case 



of degenerate density operators. A kinematic approach to 
define GP in mixed states undergoing nonunitary evolu- 
tion, generalizing the results of the above two works, has 
recently been proposed by Tong et al. [10 . Wang et al. 
[TTlfT2] defined a GP based on a mapping connecting den- 
sity matrices representing an open quantum system, with 
a nonunit vector ray in complex projective Hilbert space, 
and applied it to study the effects of a squeezed-vacuum 
reservoir on GP. 

The geometric nature of GP provides an inherent fault 
tolerance that makes it a useful resource for use in devices 
such as a quantum computer [T3| . There have been propos- 
als to observe GP in a Bose-Einstein-Josephson junction 
[14j and in a superconducting nanostructure |15j . and of 
using it to control the evolution of the quantum state [TB] . 
However, in these situations the effect of the environment 
is never negligible [17] . Also in the context of quantum 
computation, the qubits are never isolated but under some 
environmental influence. Hence it is imperative to study 
GP in the context of Open Quantum Systems. An impor- 
tant step in this direction was taken by Whitney et al. 
[18] . who carried out an analysis of the Berry phase in a 
dissipative environment |19| . Rezakhani and Zanardi |20| 
and Lombardo and Villar [21] have also carried out an 
open system analysis of GP, where they were concerned, 
amongst other things, with the interplay between decoher- 
ence and GP brought about by thermal effects from the 
environment. Sarandy and Lidar (22] have introduced a 
self-consistent framework for the analysis of Abelian and 
non-Abelian geometric phases for open quantum systems 
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undergoing cyclic adiabatic evolution. The GP acquired 
by open bipartite systems has recently been studied by Yi 
et al. [23] using the quantum trajectory approach. 

In this paper we make use of the method of Tong et al. 
[To] to study the GP of a qubit (a two-level quantum sys- 
tem) interacting with different kinds of system-bath (en- 
vironment) interactions, one in which there is no energy 
exchange between the system and its environment, i.e., 
a quantum non-demolition (QND) interaction and one in 
which dissipation takes place [24'.' 2 5]. Throughout, we as- 
sume the bath to start in a squeezed thermal initial state, 
i.e., we deal with a squeezed thermal bath. The physical 
significance of squeezed thermal bath is that the decay rate 
of quantum coherences in phase-sensitive (i.e., squeezed) 
baths can be significantly modified compared to the decay 
rate in ordinary (phase-insensitive) thermal baths |26ji27, 
[55] . A method to generate GP by making use of a squeezed 
vacuum bath has recently been proposed by CaroUo et al. 
[29]. 

The open system effects studied below can be given an 
operator-sum or Kraus representation [30] . In this repre- 
sentation, a superoperator £ due to environmental inter- 
action, acting on the state of the system is given by 



P 



k 



ek\U{p®\h){h\)U^\eu)=Y,E,pE], 



(1) 

where U is the unitary operator representing the free evo- 
lution of the system, reservoir, as well as the interaction 
between the two, {|/o)} is the environment's initial state, 
and {|efc)} is a basis for the environment. The environ- 
ment and the system are assumed to start in a separa- 
ble state. In the above equation, Ej = (efc|C/|/o) are the 
Kraus operators, which satisfy the completeness condition 
Y,^ E]Ej = I. The operator sum representation is not 
unique. Every (infinitely many) possible choice of tracing 
basis {|efe)} in Eq. ([T]) yields a different, but equivalent 
and unitarily related, set of Kraus operators. It can be 
shown that any transformation that can be cast in the 
form ll]) is a completely positive (CP) map f31]. 

From the viewpoint of quantum communication, these 
open quantum system effects correspond to noisy quan- 
tum channels, and are recast in the Kraus representation. 
We find that some of them may be interpreted in terms 
of familiar noisy quantum channels. This abstraction will 
enable us to connect noisy channels directly to their effect 
on GP, bypassing system-specific details. Visualizing the 
effect of these channels on GP in a Bloch vector picture of 
these open system effects helps to interpret our GP results 
in a simple fashion. 

The structure of the paper is as follows. In Section [2l 
we briefly discuss QND open quantum systems and collect 
some formulas which would be of use later. In Section [3l 
we study the GP of a two-level system in QND interac- 
tion with its bath. Here we consider two different kinds of 
baths. In Section [3Tl a bath of harmonic oscillators is con- 
sidered, and we also briefly touch upon a bath of two-level 
systems. In Section [X^ we point out that the GP results 
obtained in this section are generic for any purely dephas- 



ing channel. In Section (H we study the GP of a two- level 
system in a dissipative bath. Section [4 . II considers the sys- 
tem interacting with a bath of harmonic oscillators in the 
weak Born-Markov, rotating- wave approximation (RWA). 
In Section 221 "we point out that the GP results obtained 
in this section are generic for any squeezed generalized am- 
plitude damping channel [32 , of which the familiar gen- 
eralized amplitude damping channel [31j is a special case. 
We make our conclusions in Section [5] 



2 QND open quantum systems - A 
recapitulation 

To illustrate the concept of QND open quantum systems 
we use the percept of a system interacting with a bath of 
harmonic oscillators. Such a model, for a two-level atom, 
has been studied [551IMII55] in the context of influence of 
decoherence in quantum computation. We will consider 
the following Hamiltonian which models the interaction 
of a system with its environment, modelled as a bath of 
harmonic oscillators, via a QND type of coupling ^28j 

El = Hs + Hr + HsB^ 

= Hs + Yl ^^kblbk + Hs^ gk{bk + bl) 



HIE 



Tiojk 



(2) 



Here Hs, Hj^ and i?si? stand for the Hamiltonians of 
the system (S), reservoir (i?) and system- reservoir (S- 
R) interaction, respectively. The last term on the right- 
hand side of Eq. (1) is a renormalization inducing 'counter 
term'. Since [Hs,Hsr] = 0, (1) is of QND type. Here Hs 
is a generic system Hamiltonian which we will use in the 
subsequent sections to model different physical situations. 
The system plus reservoir complex is closed obeying a uni- 
tary evolution given by 



Pit) 



(3) 



where p(0) = ^^^(0)^^(0), i.e., we assume separable initial 
conditions. Here we assume the reservoir to be initially in 
a squeezed thermal state, i.e., a squeezed thermal bath, 
with an initial density matrix /9_r(0) given by 



PR 



(0) ^ Sir,<P)pt,,S\r,<P), 



exp 



where pth = Ilfc [l 
sity matrix of the thermal bath, and 

S{rk,<Pk) = exp 



(^-Phujkblbi 



(4) 

is the den- 




2 



is the squeezing operator with , being the squeezing 
parameters [36j . In an open system analysis we are inter- 
ested in the reduced dynamics of the system of interest 
S which is obtained by tracing over the bath degrees of 



Subhashish Banerjee, R. Srikanth: Geometric phase of a qubit interacting with a squeezed-thermal bath 



3 



freedom. Using Eqs. ^ and ^ and tracing over the bath 
we obtain the reduced density matrix for S, in the system 
eigenbasis, as [28] 



Here 

and 

lit)-- 



E 



sin(wfet), 



(5) 
(6) 



9l 



iE^coth(^)|(e- 



l)cosh(rfc) 



1) sinh(rfc)e 



(7) 



For the case of an Ohmic bath with spectral density I{uj) = 
^a;e~"/"=, where 70 and uJc are two bath parameters, r]{t) 
and 7(4) have been evaluated in [28], where we have for 
simplicity taken the squeezed bath parameters as 

cosh(2r(a;)) — cosh(2r'), sinh (2r(ijj)) = sinh(2r), 
<l^(uj) — auj, 

with a being a constant depending upon the squeezed 
bath. We will make use of Eqs. (O and ([7]) in the sub- 
sequent analysis (cf. Ref. [28j for details). Note that the 
results pertaining to a thermal bath can be obtained from 
the above equations by setting the squeezing parameters 
r and (i.e., a) to zero. 



3 GP of two- level system in QND interaction 
with bath 

In this section we study the GP of a two-level system in 
QND interaction with its environment (bath) . We consider 
two classes of baths, one being the commonly used bath of 
harmonic oscillators [H], and the other being a localized 
bath of two-level systems. 

3.1 Bath of harmonic oscillators 

The total Hamiltonian of the S + R complex has the 
same form as in Eq. ([2]) with the system Hamiltonian 
Hs = ^cr^, where is the usual Pauli matrix. We will be 
interested in obtaining the reduced dynamics of the sys- 
tem. This is done by studying the reduced density matrix 
of the system whose structure in the system eigenbasis is 
as in Eq. ((Sj) . For the system described by Hs an appropri- 
ate eigenbasis is given by the Wigner-Dicke states [5711^ 
[39] I j, m) , which are the simultaneous eigenstates of the 
angular momentum operators and Jz, and we have 
Hs\j,m) = hujm\j,m) = Ej^rn\ji'rri) . Here —j<m< j. 
For the two-level system considered here, j — ^ and hence 
m = Using this basis in Eq. ^ we obtain the re- 

duced density matrix of the system as 



Pjm^jn (^) 



^ — iuj{in—ri)t^i{Tiuj) {ra —n )ri{t) 



(8) 



It follows from Eq. ([8]) that the diagonal elements of the 
reduced density matrix signifying the population remain 
unaffected by the environment whereas the off-diagonal 
elements decay. This is a feature of the QND nature of 
the system-environment coupling. Initially we choose the 
system to be in the state 

1^(0)) =cos(|)|l)+e^*« sin (^l)lO). (9) 

Using this we can write Eq. ^ as 



Pjojoit) 


— cos^( 


2 ' 




PkM^) 


= lsin 






Piuoit) 


^Isin 


Jgg)gi('^t+0o)g-(?iw)^7(t) 




Pkiii^) 


= sin2( 


2 ^ 


(10) 



We will make use of Eq. pU]) to obtain the GP of the above 
open system using the prescription of Tong et al. [lOj 



N 



$GP = arg ^ VAfc(0)Afc(T)(tf'fc(0)|tf'fe(T)) : 
\fc=i 
-/^"<it(*,(t)|*'fc(t))' 



e Jo 



(11) 



Hereafter we will consider for GP a quasi-cyclic path where 
time (t) varies from to r = 27r/ti;, w being the system 
frequency. In the above equation the overhead dot refers 
to derivative with respect to time and Afe(r), !Z'fe(r) refer 
to the eigenvalues and the corresponding eigenvectors, re- 
spectively, of the reduced density matrix given here by Eq. 
([To]) . The eigenvalues of Eq. (fTO]l are 



1 



A±(t) = -[l-Hcos(0o)e±(O] 



(12) 



where e±{t) = ±y^l + tan2(6lo)e-2(M=^7(«). Since 7(0 = 
for t = 0, we can see from the above equations that 
A+(0) = 1 and A_(0) = 0. From the structure of the 
Eq. ([TT]) we see that only the eigenvalue A+ and its corre- 
sponding eigenvector need be considered for the GP. 
This normalized eigenvector is found to be 

|f+(t)) =sin('|] |l)+e'M+*°)cos('|) |0), (13) 



where sin {Ot/2) 



e+ + l 



It can be seen that for < = 0, 



sin (y) cos (^) and co s (-^ ) sin (^), as expected. 
Now we make use of Eqs. ([T2J, ([13]) in Eq. ([TT]) to obtain 
GP as 



^GP = arg 



1 + cos(e'o)y 1 + tan2(6'o)e-2(M^7(r) 
X |cos(y)sm l—j +e"^^sm i—j cos(y)| 

^ g-»^/;<itcos^(%)i_ ^^^^ 
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Here '-f{t) is as given in Ref. ([2H]) for a zero temperature 
(T) bath or high T bath. It can be easily seen from Eq. 
(fT4|) that if we set the influence of the environment, en- 
capsulated here by the expression j{t), to zero, we obtain 
for T = ^, -Pop = -^/^ = -tt{1 - cos(6'o)), where n is 
solid angle subtended by the tip of the Bloch vector on the 
Bloch sphere, which is the standard result for the unitary 
evolution of an intial pure state. More generally, unitary 
evolution of mixed states also has a simple relation to the 
solid angle, given by 

<Pgp = - tan-i (^L tan ^ j , (15) 

where L is the length of the Bloch vector 

The effect of temperature and squeezing on GP is broug' 
out by Figs. [Hand [21 From Figs.[HA) and (B), we see, re- 
spectively, that increasing the temperature and squeezing 
induce a departure from unitary behavior by suppressing 
GP, except at polar angles 6*0 = 0, 7r/2 of the Bloch sphere. 
It can be shown that, similarly, increase in the S-R cou- 
phng strength, modelled by 70, also tends to suppress GP. 
(Throughout this article, the Figures use uj = 1. Further, 
Figures in this Section use Wc = 40^.) The suppresive 
influence of temperature on GP is also seen in Figs. [21 
where temperature is varied for fixed 9o and squeezing. A 
similar suppresive influence of squeezing on GP is brought 
out by comparing Figs.[2I^A) and[2tB). These observations 
are easily interpreted in the Bloch vector picture, as we 
discuss later in this section. 

Another interesting case is that of qubit subjected to 
a bath of two-level systems, studied by Shao and collabo- 
rators in the context of QND systems [4^, and quantum 
computation [41]. It has also been used to model a nano- 
magnet coupled to nuclear and paramagnetic spins [32] ■ It 
can be shown 43J that this case is mathematically similar 
to that of QND interaction with a vacuum bath of har- 
monic oscillators for weak S-R coupling, and hence the 
dependence of GP on 6q and 70 is similar to the analo- 
gous case discussed above. 




Fig. 1. GP (Eq. (O) as a function of 60 (in radians) for 
different temperatures and squeezing at 70 — 0.0025. In both 
plots, unitary evolution is depicted by the large-dashed curve. 
(A) GP at r = a = 0.0; the dot-dashed, small-dashed and 
solid curves correspond, respectively, to temperatures 50, 100, 
300. (B) GP at r = 100 and a = 0; the dot-dashed, small- 
dashed and solid curves correspond, respectively, to squeezing 
parameter r = 0, 0.4, 0.6. For QND interactions, in the region 
tt/2 < Oq < n, the pattern is symmetric but sign reversed. 
Observe that, as is true for all QND cases, GP vanishes at 
9o = 0. This can be attributed to the fact that the qubit 's 
evolution sweeps no solid angle in this case. Here, as in all 
other Figures, we take oj = 1, and for all Figures in this Section, 
LUc = 40a;. 



3.2 Evolution of GP in a phase damping channel 

While the results derived above are for QND S-R inter- 
actions with two types of baths, they are quite general, 
and in fact apply to any open system effect that can be 
characterized as a phase damping channel [31j . This is 
a uniquely non-classical quantum mechanical noise pro- 
cess, describing the loss of quantum information without 
the loss of energy. This system can be represented by the 
Kraus operator elements 



1 

e 







;/3(t 



(16) 

where (3(t) encodes the free evolution of the system and 
X{t) the effect of the environment. It is not difficult to see 
that the QND interactions we have considered realize a 
phase damping channel. 



In the case of QND interaction with a bath of harmonic 
oscillators (Sec. 13. ip . it is straightforward to verify that 
with the identification 

X{t) = 1 - exp [-2{nwfj{t)] ; P{t) = cot. (17) 

the operators ([16]) acting on the state ([9]) reproduce the 
evolution Eq. ([T0|) by means of the map Eq. ([1]) . Similarly, 
the effect of QND interaction with a bath of two level sys- 
tems can also be represented as phase damping channel 
[33]. Our result is in agreement with that of Ref. |llj . 
where GP is shown to depend on the dephasing parame- 
ter, introduced phenomenologically. Our result is obtained 
from a microscopic model, governed by Eqs. ([21)-((31), that 
takes into consideration the interaction of a qubit with 
a squeezed thermal bath, the resulting dynamics being 
shown above to be equivalent to a phase damping chan- 
nel. 
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Fig. 2. GP (in radians) as a function of temperature (T, in 
units where h = ks = 1) for QND interaction with a bath 
of harmonic oscillators (Eq. (|14|) '). (A) with 70 — 0.005 and 
vanishing squeezing. The solid, dashed and larger-dashed lines 
correspond to ^0 ~ vr/8,37r/16 and 7r/4. (B) Same as Figure 
(A), except that here squeezing is non- vanishing, with r — 0.7 
and a — 0.1. 



In the case of QND interaction, any initial state not lo- 
cated on the cr3-axis tends to inspiral towards it, its trajec- 
tory remaining coplanar on the x-y plane. Consequently, 
the entire Bloch sphere shrinks into a prolate spheroid, 
with its axis of symmetry given by the (T3 axis. The extent 
of inspiral depends upon the parameter X{t); the greater 
is X(t), the more is the inspiral. Greater squeezing and 
higher temperature accentuate this shrinking. 

Guided qualitatively by the relation Eq. we may 
interpret GP as directly dependent on the Bloch vector 
length L{t), and the sohd angle (]7) subtended at the cen- 
ter of the Bloch sphere during a cycle in parameter space. 
Increasing T, 70 or squeezing results in a larger degree of 
inspiral causing a reduction of both L and i7, and hence 
greater suppression of GP relative to the case of unitary 
evolution. 

In Figs.[lJA) and (B), we noted that the GP remains 
invariant at polar angles do — and Oq — 7r/2. In the 
case 6*0 = 0, the Bloch vector remains a constant (0, 0, 1) 
throughout the evolution and hence accumulates no GP. 
In the case 6*0 = 7r/2, note that SI = 2tt. From Eq. (fT5|) . we 
see that irrespective of the length of the Bloch vector, GP 
should remain the same, i.e., — tt. This suggests that in the 
general nonunitary case, when the Bloch vector rotates on 



the equitorial plane, GP is unaffected by whether or not 
there is an inspiral of the Bloch vector. 

The fall of GP as a function of T (Figs.^A) and© can 
be attributed to the fact that as T increases the tip of the 
Bloch vector inspirals more rapidly towards the (73 axis, 
and thus sweeps less GP. Squeezing has the same effect as 
temperature, of contracting the Bloch sphere along the CT3 
axis, leading to further suppression of GP (Figs. [T](B) and 
HB)). 



4 GP of two- level system in non-QND 
interaction with bath 

In this section we study the GP of a two-level system in a 
non-QND interaction with its bath which we take as one 
composed of harmonic oscillators. We consider the case of 
the system interacting with a bath which is initially in a 
squeezed thermal state, in the weak coupling Born-Markov 
RWA. 



4.1 System interacting with bath in the weak 
Born- Markov RWA 

Now we take up the case of a two-level system interacting 
with a squeezed thermal bath in the weak Born-Markov, 
rotating wave approximation. This kind of system-reservoir 
(S—R) interaction is consonant with the realization that 
in order to be able to observe GP, one should be in a 
regime where decoherence is not predominant p!51[^ . The 
system Hamiltonian is Hs and it interacts with the bath of 
harmonic oscillators via the atomic dipole operator which 
in the interaction picture is given as 



D{t) = da^ 



J* _ i\ 



(18) 



where d is the transition matrix elements of the dipole 
operator. The evolution of the reduced density matrix op- 
erator of the system S in the interaction picture has the 
following form [44ll45] 



-p^(0=7o(iV + l) 

+ loN (^a+p'{t)a_ - ^<J-(J+pHt) - ^p-{t)<j_a+ 
~ "foMa+p%t)a+ - 7oA/V_p"(t)a_. 



(19) 



Here 70 is the spontaneous emission rate given by 70 = 
Auj^\d\'^ /3tic^ , and a+, ct_ are the standard raising and 
lowering operators, respectively given by 

a+ = |1) (0| = ^ (^1 + ^^2) ; a_ - |0>(1H ^ {a, - za^) . 

(20) 
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Eq. 
as 



\l9\i may be expressed in a manifestly Lindblad form 



the purpose of GP we need only the eigenvalue A+ {t) , and 
its corresponding normalized eigenvector is given as 



), (21) |f+W)=sin(^|)|l) 



e'(x(t)+c^t) cos 



i ) 10)' 



(27) 



= (7o7Vth/2)i/^i?t 
This observation 



where i?i = (7o(iVth + 1)I2Y'^R. i?2 
and R = a- cosh(r) + e**(T+ sinh(r) 
guarantees that the evolution of the density operator can 
be given a Kraus or operator-sum representation [3T], a 
point we return to later below. If T = 0, then i?2 van- 
ishes, and a single Lindblad operator suffices to describe 

Eq. (ng). 

In the above equation we use the nomenclature |1) for 
the upper state and |0) for the lower state and cti, (72, cts 
are the standard Pauli matrices. In Eq. (|19p 



N = 7Vth(cosh^(r) + sinh^(r)) + sinh^(r), 
M = -isinh(2r)e^*(2iVth-f 1), 



1 



(22) 



Here iVth is the Planck distribution giving the number of 
thermal photons at the frequency m and r, <P are squeezing 
parameters. The analogous case of a thermal bath without 
squeezing can be obtained from the above expressions by 
setting these squeezing parameters to zero. We solve the 
Eq. (|19p using the Bloch vector formalism to obtain the 
reduced density matrix of the system in the Schrodinger 
picture as [43] 



2 



5*e'"* ^1-^) 



(23) 



where, 



A^(a3(i))-e-^''^''^+^)*(a3(0)) 



(2N+1) 



B 



-7o(2Af+l)t 



-^(27V+l+a)t 



(24) 



-I- sinh( 



(p^o"^ - i ) e 
^°"^e'*-^(2^+i)*(a+(0)). 



a_(0)) 



(25) 



Here a = sinh(2r)(2iVth + 1). Making use of Eq. Eq. 
(|25|) can be written as B = Re~'^^. The explicit expres- 
sions for R and x may be found in Ref . 43J . For the deter- 
mination of GP we need the eigenvalues and eigenvectors 
of the Eq. ([23]) . The eigenvalues are 



A±(t) = ^(l + e±) 



(26) 



where e± — ±V + AR"^. As can be seen from the above 
expressions, at t = 0, A+(0) = 1 and A_(0) = 0, hence for 



27?. 



It can be 



where sin((?f/2) , 

seen that for t = 0, x(0) = 0o, sin ) = y^i+I^M ^ 

cos(f) and cos(f) = ^ll^Im. = sin(|i), as ex- 
pected. Now we make use of Eqs. ([Ml), dSZ]) in Eq. HI]) to 
obtain GP as 

«?GP = arg ' 



{-(l + yi^(7yT4i?^)}^ 



X < cos I — I sm 



X e 



i ^ dt(x(t)+u,)cos\^.) 



(28) 



It can be easily seen from the Eq. ([28| that if we set the 
influence of the environment, encapsulated here by the 
terms 70, a and to zero, we obtain for t = = 
— 7r(l — cos(0o))i as expected, which is the standard result 
for the unitary evolution of an intial pure state [8l|9] . Thus 
we see that though the Eqs. (fT4]) . (fSS]) represent the GP 
of a two-level system interacting with different kinds of S- 
R interactions, when the environmental effects are set to 
zero they yield identical results. This is a nice consistency 
check for these expressions. 

As expected, increasing the temperature, S—R cou- 
pling strength or squeezing induces a departure of GP 
from unitary behavior. However the interpretation is less 
straightforward than in the QND case. Further, introduc- 
tion of squeezing complicates this pattern by disrupting 
the monotonicity of the GP plots, as evident from the 
'humps' seen for example in the Fig. EjjB), in comparison 
with those in Fig.[3]^A). 

In all cases, we find that GP vanishes at 6q — tt, i.e., 
for a system that starts in the south pole of the Bloch 
sphere. On the other hand, for sufficiently small 70, we 
find from Figs.[3][A) and[3]^B) that GP may vanish also in 
the case 9^ = 0. These observations may be interpreted in 
the Bloch vector picture, and are discussed in Section [521 

In contrast to the situation in a purely dephasing sys- 
tem, GP in a dissipative system is rather complicated, and 
less amenable to interpretation. The dependence of GP on 
temperature is depicted in Figs. [5] and [3 The expected 
pattern of GP falling asymptotically with temperature is 
seen. Our results parallel those obtained in Refs. [20ll46] 
for the case of zero squeezing (Figs. [U^A) and[HA)), and 
extend them to the case of a squeezed thermal environ- 
ment. We note that the effect of squeezing is to make GP 
vary more slowly with temperature, by broadening the 
peak and fattening the tails of the plots. This counterac- 
tive behavior of squeezing on the influence of temperature 
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Fig. 3. GP as a function of 60 (in radians) for different val- 
ues of 70 and squeezing in the Born-Markov approximation 
(Eq. (|28|) ). The discontinuity in GP after tt is due to the con- 
vention that an angle in the third quandrant is treated as 
negative. (A) T = 0. The large-dashed curve is the unitary 
case (70 = 0) . The dot-dashed (small-dashed) curve represents 
70 = 0.1 (70 = 0.3). The solid curve represents 70 = 0.6. 
The stationary state, for which GP vanishes, corresponds to 
do = n (i.e., |0)), to which all states in the Bloch sphere are 
asymptotically driven. Thus, a qubit started in this state re- 
mains stationary and acquires no GP. (B) Same as Figure (A), 
except that squeezing r = 0.4, <1> — 7r/4. 



on GP for the case of a dissipative system is interesting, 
and would be of use in practial implementation of geomet- 
ric phase gates. This effect can be understood by visualiz- 
ing the effects of squeezing and temperature on the Bloch 
sphere, a point we return to in Section l472l 



4.2 Evolution of GP in a squeezed generalized 
amplitude damping channel 

While the results derived in this section pertain to a dissi- 
pative S-R interaction in the Born-Markov RWA, they are 
quite general, and are applicable to any open system effect 
that can be characterized as a squeezed generalized ampli- 
tude damping channel [32j . Amplitude damping channels 
capture the idea of energy dissipation from a system, for 
example, in the spontaneous emission of a photon, or when 
a spin system at high temperature approaches equilibrium 
with its environment. A simple model of an amplitude 
damping channel is the scattering of a photon via a beam- 
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Fig. 4. GP (in radians) vs temperature (T, in units where 
h = kB = l) from Eq. (HHJ. Here w = 1.0, 6*0 = 7r/2, the large- 
dashed, dot-dashed, small-dashed and solid curves, represent, 
respectively, 70 = 0.005, 0.01, 0.03 and 0.05. (A) squeezing 
is set to zero; (B) squeezing non-vanishing, with r = 0.4 and 
= 0. 



splitter. One of the output modes is the environment, 
which is traced out. The unitary transformation at the 
beam-splitter is given by i? = exp [9{a^b ~ a6^)] , where 
a, h and a^, 6^ are the annihilation and creation operators 
for photons in the two modes. The generalized amplitude 
damping channel, with T > and with zero squeezing, 
extends the amplitude damping channel to finite tem- 
perature [31j . A very general CP map generated by Eq. 
(|T9| has been recently obtained by us [32], and could be 
appropriately called the squeezed generalized amplitude 
damping channel. This extends the generalized amplitude 
damping channel by allowing for finite bath squeezing. It 
is characterized by the Kraus operators [32] 



-Bo = 

E2 = 

E3 = 



1 




v/l-M(t)0 

VI - 



VT^We-^* 



(29) 
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Fig. 5. GP vs temperature (T, in units where ft = fes = 1) 
from Eq. (I28j. Here tj = 1.0, 9o = 7r/2 + 7r/4. The curves 
represent 70 = 0.005, 0.01, 0.03 and 0.05 as in Fig. H (A) 
squeezing is set to zero; (B) squeezing non- vanishing, with r = 
0.4 and = 0. 



With some algebraic manipulation, it can be verified that 
with the identification 



N 



P2(2iV + l 



-7o(2Ar+l)t\ 



2N+1 sinh^(7oaV2) 



, , _i_expf-^(2iV+l)t 
2p2iV sinh(7o(2Ar + l)t/2) 2^ 

(t) = - (1 - P2 W) + y{t)\ - e-^°(2iv+i)*) , (30) 



where N is as in Eq. ((22| . the operators (p9|) acting on 
the state © reproduce the evolution (|23p . by means of 
the map Eq. ([1]), provided p2 = 1 — Pi, satisfies 



1 



P2 = 



41? 



(31) 



X [A^S + + AiB^ - C - B{l + C) - D) 

- (1 + B)!) - C(B + £> - 1) 

± 2 (i:>(B - + (yl - 1)C + D) 

X (A-AB + (B-l)C + i:'))^/2 



where 



2iV+l sinh^(7oaV2) 
2N sinh(7o(27V + l)t/2) 



exp(-7o(2iV+l)V2), 



N 

^ ^ 2ivTT^^ " cxp(-7o(2iV + l)t)), 

C = A + B + cxp(-7o(2iV + l)t), 

D = cosh2(7oai/2) exp(-7o(27V + l)t). 



(32) 



As the interaction in the Born-Markov RWA realizes 
a squeezed generalized amplitude damping channel [32], 
the various qualitative features of GP seen under a dissi- 
pative interaction (for example, the relatively complicated 
dependence of GP on and on evolution time) carry over 
to any squeezed generalized amplitude damping channel. 
If squeezing parameter r is set to zero, it can be seen from 
above that Eq. ([^H) reduces to a generalized amplitude 
damping channel, with iy{t) = a(t), fi{t) = and pi and 
P2 being time-independent. If further T = 0, it can be 
seen from above that p2 = 0, reducing Eq. (^5)) to two 
Kraus operators, corresponding to an amplitude damping 
channel. 

Refs. [TT] and [H] consider GP evolving under an am- 
plitude damping channel and a squeezed amplitude damp- 
ing channel, respectively. These are subsumed under the 
squeezed generalized amplitude damping channel consid- 
ered above. This channel is contractive, in that the system 
is seen to evolve towards a fixed asymptotic point in the 
Bloch sphere, which in general is not a pure state, but the 
mixture 

Pasymp ~ I Q q I ' ^ ' 

where q= (A^+l)/(2iV+l). If T = r = 0, then g = 1, and 
the asymptotic state is the pure state |0). Physically this 
can be understood as a system going to its ground state 
by equilibriating with a vacuum bath. This can have a 
practical application in quantum computation in the form 
of a quantum deleter [37]. At T ^ 00, p ^ 1/2, and the 
system tends to a maximally mixed state, thereby realizing 
a fully depolarizing channel 31j. 

As in the case of the QND interaction, abstracting the 
effect of dissipative interaction into the Kraus represen- 
tation allows us to subsume all the details of the system 
into a limited number of channel parameters pi (i), <P, a(t), 
fi(t) and 1^(1). Any other dissipative system that can be 
described by a Lindblad-type master equation Eq. 
will show a similar pattern in behavior. 

To develop physical insight into the solution, we trans- 
form to the interaction picture, and for simplicity, set the 
squeezing parameters to zero. Then, the action of the oper- 
ators (|29p . [which now represents a generalized amplitude 
channel] on an arbitrary qubit state is given in the Bloch 
vector representation by 



('T(i)) 



((ai(o))yT^A(t), {a^m^T^m, 

A(t)(l-2p) + (a3(0))(l-A(t))), (34) 



where p = {Nti, + l)/(27Vth + 1) and X{t = 00) = 1. Thus, 
the Bloch sphere contracts towards the asymptotic mixed 
state (0,0, 1 — 2p) (Fig. [6lJA)), characteristic of a gener- 
alized amplitude damping channel, with T > and no 
squeezing. If T = case, then p = 1, and the asymptotic 
state (0,0,-1) is pure. 
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(A) 





([24| that for sufficiently large 70, (o'sit)) changes sign at 
ti EE log(2[iV + l])/(7o[27V + 1]). Further, we note that R 
vanishes for Oq — (as well as 9q = tt) . 

It is convenient to recast Eq. ([28|) in the expanded form 

<I>GP = tan-1 [(sin(x(T) - x(0) + 2n) sin(0o/2) cos(0,/2)) 
- {cos(x(t) - x(0) + 2tt) sin(0o/2) cos(0,/2) 
+ cos(6lo/2)sin(6i^/2)}] 

(35) 



Fig. 6. Shrinking of the full Bloch sphere into an oblate 
spheroid under evolution given by a Born-Markov type of dis- 
sipative interaction with 70 = 0.6 and temperature T — 5.0. 
In (B), the x-y axes are interchanged for convenience. (A) 
r = <P ^ 0, t ^ 0.15; (B) r = OA,<P = 1.5, t = 0.15. Finite 
(l> is responsible for the tilt. 



The Bloch vector picture allows us to interpret the re- 
sults of Section HTT] Eqs. show that the Bloch vector 
for the states corresponding to = 0, tt move only along 
the 2;-axis of the Bloch sphere for zero as well as finite T. 
For the case 9q = tt and zero T, the Bloch vector remains 
stationary at (0, 0, —1), and hence GP vanishes. In the fi- 
nite T case, GP still vanishes, because the Bloch vector has 
the form (0, 0, —L{t)), where the Bloch vector length L{t) 
shrinks from 1 towards an interaction-dependent asymp- 
totic value, which is zero for infinite temperature or finite 
otherwise. Since the Bloch vector shrinks strictly along 
its length, and thus subtends no finite angle at the cen- 
ter of the sphere, we find that GP vanishes at 6*0 — tt, as 
expected (cf. Figs.E]). 

On the other hand, even though the Bloch vector shrinks 
similarly along its length in the case 6*0 = 0, we find that 
GP is non-vanishing in certain cases, in fact, in precisely 
those cases where the tip of the Bloch vector crosses the 
center of the Bloch sphere moving along the ds-axis. That 
is, they correspond to the situation where (0-3 (t)) changes 
sign from positive to negative during the period of one 
cycle. In these cases, the dependence of GP on the Bloch 
vector is too involved for us to interpret in terms of L 
and the angle subtended by the Bloch vector, for some 
qualitative insight. Nevertheless this feature may be for- 
mally understood as follows. It can be observed from Eq. 



dtixit) + u) cos' ^ 



It can be seen that for the case = cos{9t/2) — 1 and, 
in particular, cos{9r/2) — 1. Substituting these values in 
Eq. p5|) , it is seen that GP vanishes because the two terms 
in the RHS of Eq. ((35|) cancel each other. Next consider 
the case where = but where 70 is sufficiently weak 
that T < ti, i.e., {(Tsit)) does not change sign during one 
cycle. In this case, from above it is seen that cos(0(/2) = 0, 
and, in particular, cos(0r/2) — 0, and thus the terms in the 
RHS of Eq. ((35| vanish identically. But in the case of 6*0 — 
where t > ti (70 being relatively stronger), cos(0t/2) = 
initially in the time interval [0, ti] , and then switches to 1 
in the interval (ti, r]. In particular, cos(6'i-/2) = 1. Observe 
that if cos(0t/2) — 1 throughout the interval [0, t], the two 
terms in the RHS cancel each other. It follows that GP is 
non- vanishing because of an excess contributed by the first 
term, in the interval [0, ii]. 

Contraction produced by an increase in temperature 
tends to be less pronounced in the presence (than in the 
absence) of squeezing (Figs. [5]). This is reflected in the 
slower variation of GP with respect to temperature, seen 
in Figs.HjB) andEl^B) in relation to Figs.[l^A) andElJA), 
respectively. As observed in Figs. H] and [HI GP falls as 
a function of T, for sufficiently large T. This may quite 
generally be attributed to the reduction in L and f2 caused 
by the contraction of Bloch vector as a result of interaction 
with the environment. The tilt of the contracted Bloch 
sphere in Fig. \ElB) is due to finite <!>. 



5 Conclusions 

We have studied the combined influence of squeezing and 
temperature on the GP for a qubit interacting with a bath 
both in a non-dissipative as well as in a dissipative man- 
ner. In the former case, squeezing has a similar debilitat- 
ing effect as temperature on GP. In contrast, in the latter 
case, squeezing can counteract the effect of temperature 
in some regimes. This makes squeezing potentially helpful 
for geometric quantum information processing and geo- 
metric computation. In particular, in the context of using 
engineered (e.g., squeezed) reservoirs to generate GP 
it would be helpful to consider the effect of squeezing to- 
gether with thermal effects p0 l l21 ) . 

In the non-dissipative (QND) case, we analyzed a num- 
ber of open system models using two types of bath: the 
usual one of harmonic oscillators, and that of two-level 
systems. It was shown that for the case of weak S—R cou- 
pling, the two kinds of baths can be mapped onto each 



10 



Subhashish Banerjee, R. Srikanth; Geometric phase of a qubit interacting with a squeezed-thermal bath 



other. GP was studied as a function of the initial polar 
angle 0o of the Bloch sphere, temperature and squeez- 
ing (arising from the squeezed thermal bath) . In the QND 
case, it was seen that increasing 70 , temperature or squeez- 
ing tends to cause a similar departure from unitary behav- 
ior by suppressing GP. 

However, in the dissipative case (with the environment 
modelled as a squeezed thermal bath in the weak Born- 
Markov RWA), we found that the dependence of GP on 
^0, temperature and squeezing shows a greater complex- 
ity. Here, an interesting feature due to squeezing is that 
it can disrupt, over an interval, the otherwise monotonic 
behavior of GP as a function of (the humps seen in 
Figure[3KB)). More pronouncedly, the counteractive effect 
of squeezing on temperature is brought out by a compari- 
son of FiguresHJA) withHKB), andE^A) withE^B). Also, 
its effect on the Bloch sphere is to shrink it to an oblate 
spheroid, in contrast to a QND interaction, which pro- 
duces a prolate spheroid. Thus, an interesting feature that 
emerges from our work is the contrast in the interplay be- 
tween squeezing and thermal effects in non-dissipative and 
dissipative interactions. By interpreting the open quan- 
tum effects as noisy channels, we make the connection 
between geometric phase and quantum noise processes fa- 
miliar from quantum information theory. 

An added feature of our work is that we make a con- 
nection between the studied open system models and the 
phase damping and the newly introduced squeezed gen- 
eralized amplitude damping [32j channels, noise processes 
which are important from a quantum information theory 
perspective. In particular, we give a detailed microscopic 
basis for these noisy channels. This allows us to study the 
effects of the formal noise processes on GP. 
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